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ABSTRACT 

It is proposed that the gamma ray burst photons near the peak of the spectrum at several 
hundred KeV are produced on very compact scales, where photon production is limited by 
blackbody effects and/or the requirement of energetic quanta (E > 2m e c 2 ) for efficient further 
production. The fast variation of order milliseconds in the time profile is then a natural 
expectation, given the other observed GRB parameters. Analytic calculations are presented 
to show that the escape of non-thermal, energetic gamma rays can emerge within a second of 
the thermal photons from a gammasphere of below 10 12 cm. The minimum asymptotic bulk 
Lorentz factor in this model is found to be of order several hundred if the photosphere is of 
order 3 x 10 11 cm and greater for larger or smaller photospheric radii. It is suggested that 
prompt UHE gamma rays might provide a new constraint on the asymptotic Lorentz factor of 
the outflow. 



1. Introduction 

Gamma ray burst spectra typically peak at several hundred KeV, and frequently (though not always) 
have an extended non-thermal tail that contains a significant fraction of the total energy. Although the 
high energy side of the peak is hard to characterize with reliable generalizations, classical GRB's seem to 
typically peak above 100 KeV; a significantly softer peak is a sign that the GRB is in a separate class of 
"soft repeaters". Though this has been suggested to be a selection effect, we believe it to be real. The 
total luminosity in soft X-rays is typically below the luminosity above 100 KeV by an order of magnitude 
or more. When one considers that, realistically, an emitting surface will have some dispersion in its local 
temperature (e.g. a hot spot is likely to have cool edges, relativistic beaming from material moving directly 
at the observer is likely to be mixed with contributions from fluid with a sideways component to its motion 
etc.) the X-ray paucity seems all the more clean and significant. 

We believe that the characteristic peak energy of about 200 KeV is not merely a chance value or 
selection effect, but is rather telling us something important about GRB's. In this paper we construct a 
model that, by design, yields this characteristic peak energy naturally. The basic assertion is that the 
gamma rays near the peak are emitted from a compact region of radius R not much greater than about 
10 9 rcm, where T is the bulk Lorentz factor of the expanding fireball. Earlier discussions and variations of 
this basic idea include those by Eichler (1994) and Thompson (1997). Energization of the pair plasma at 
this radius could come from internal shocks (e.g., Eichler 1994; Rees & Meszaros 1994; Sari & Piran 1997) 
but could also come from, say, collimation by surrounding baryonic material. The possibility of collimation 
may be motivated by the huge fluence from recently discovered GRB's such as GRB 990123. Interaction 
with the collimating material could produce strong internal shock, or could proceed via radiative viscosity 
with the wall of the collimating material. 

In the latter process, dissipation (i.e. creation of additional quanta) would proceed very efficiently 
when the average photon energy in the observer frame approached m e c 2 , for then large angle scattering 
by the walls of the collimating material of photons back into the jet, where they could appear in the local 
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frame to be blucshiftcd to even higher energy, would allow creation of additional pairs. Once the average 
energy per quantum in the observer frame goes far enough below m e c 2 , this dissipation mechanism would 
taper off. This is true of any dissipation mechanism that proceeds by pair creation. The point is that, if 
the clustering of GRB peaks near or just below m e c 2 is to be regarded as a signature of dissipation by pair 
production, then such dissipation would probably occur at modest T, or else the product of large T with a 
small fireball frame energy per quantum consistently yielding a given value for the observer's frame would 
seem a puzzling coincidence. If the dissipation is to take place at modest T, we argue that this probably 
indicates compact scales, before the fireball accelerates to its asymptotic value of T, which according to 
several considerations is quite large (> 300). 

At compact scales, thermalization is a strong possibility, and sets a convenient reference point for the 
analysis, but there is no implication that non-thermal effects are absent, or that the photon spectrum need 
look thermal on either side of the peak. Our view is that the non-thermal component in GRB's is not, 
insofar as can be currently established by observations, as remarkably consistent from burst to burst as is 
the spectral peak. Thus, non-thermal acceleration might proceed downstream of the photosphere, as in 
optically thin synchrotron models of bursts with only the constraint that the non-thermal component never 
dwarf the thermal component (in total photons produced). The rationale for why both components might 
be of the same order of magnitude is given in section 3, where it is shown that the energy budget of the 
burst released in the internal shocks of an unsteady flow might be distributed logarithmically over a broad 
range of radii. 

It is also possible that direct shock energization of photons (Blandford and Payne, 1982) can create 
a non-thermal photon spectrum up to 100 KeV in the frame of the fluid just below the photosphere by a 
shock that passes through the photosphere of a GRB (Eichler 1994). The photons would escape before 
they are thermalized or saturatedly Comptonized. We will argue in this paper that the material at the 
photosphere could be moving with a bulk Lorentz factor T that exceeds 300, depending on the extent 
to which non-thermal processes maintain a pair population well above the Boltzmann equilibrium value. 
Thus the high energy cutoff of the non-thermal spectrum could be above 30 MeV in the frame of the 
observer. The total photon number would have been established to within a factor of order unity well 
within the photosphere (at an optical depth of order 1/a if by bremsstrahlung) . On the other hand, the 
photon number density could change by a factor of order unity over the hydrodynamical timescale if the 
expansion, assumed to be unsteady, deviates from adiabaticity by a substantial margin. In this case, self 
absorption, which would keep the photon number steady in a static environment, fails to do so in the course 
of the expansion. As the Comptonization timescale just below the photosphere is at least as long as the 
hydrodynamical timescale, the spectrum of soft photons is not fully saturated by Comptonization. 

Since the photospheric radius in our model is considerably smaller than in some other models, which by 
design recognize the restrictions on R and T implied by the escape of energetic, non-thermal gamma rays, 
we redo the calculation of this restriction in a transparent manner and show that it is consistent with the 
other restrictions that we invoke here. These calculations show that prompt very energetic gamma rays are 
possible, even prompt TeV ones, depending on the strength of the magnetic field. Conceivably, detection 
and time resolution of very energetic gamma rays in future experiments could set a useful diagnostic of the 
field strength and bulk Lorentz factor. 

While the constraint of X-ray paucity might not be rigorous, and reconcilable with a "just-so" selection 
of parameters, we note that, historically, it was taken seriously (Imamura & Epstein 1987) before the 
cosmological GRB picture was made popular by the all sky isotropy established by BATSE observations. 
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To avoid confusion, wc define several distinct time intervals: The elapsed time in the frame of the 
fireball as the fireball evolves through radius R is called At' and is equal to R/cT. The characteristic 
timescale of variation of the emission as seen by the external observer is called 5t and is generally taken to 
be of order (in any case, at least) R/cY 2 . 

Henceforth, primed quantities denote quantities measured in the comoving frame, whereas unprimed 
quantities refer to quantities measured in the frame of the central engine. Subscripted quantities refer to 
that quantity expressed as the power of ten in cgs units denoted by the subscript. R12 for example means 
the radius R expressed in units of 10 12 cm. 

2. X-Ray Paucity 

A noteworthy feature of gamma ray burst spectra is the so-called paucity of X-rays. Although 
some BATSE spectra (though apparently not all, Preece et al. 1999) may be consistent with the 
instantaneous synchrotron radiation spectrum of a mono-energetic electron population, which would go 
as vl(ii)dv cx v 4 / 3 dv, they are often too hard to be the time integrated spectrum of such a population, 
vl(v) oc z/2 , often too hard to be the optically thin instantaneous (i.e. thin target) emission of a shock 
accelerated spectrum, vl oc , and usually too hard to be the thick target emission spectrum of shock 
accelerated electrons vl(y) oc constant. For reviews of spectral distribution see Pendleton et al. (1994) 
and Band et al. (1993). Popular fireball models often assume that the primary peak of the burst radiation 
comes from thermal electrons that have been accelerated to large enough energies to synchrotron radiate at 
several hundred KeV, but this, in our opinion, requires much fine tuning. It requires that the magnetic field 
value, the electron energy and the bulk Lorentz factor always conspire to put the peak at several hundred 
KeV; it is seldom below that for most GRB's. (It has been debated that the narrow range of vl v peaks seen 
in GRB's can be attributed to instrumental effects [Dermer, et al. 1998]. Here we are motivated to find a 
physical explanation.) Another problem with optically thin synchrotron emission is that invoking a thin 
target spectrum means low radiation efficiency, unless there is additional fine tuning. On the other hand 
allowing the electrons to radiate most of their energy in the magnetic field ("thick target emission") would 
make too much softer X-radiation. 

Invoking optically thin inverse Compton emission (IC) for the peak gamma rays would seemingly 
exacerbate the problem, because the energy of the seed photons is less likely to be monoenergetic than 
the virtual photons in the magnetic field, and even more fine tuning would be required to ensure that the 
upscattered photons always peak at ~ 300 KeV. [A saturated IC spectrum has far fewer free parameters 
(Liang 1998) and is indeed considered below, but the location of the peak still depends on the total number 
of quanta produced per unit energy. So we will consider such production first.] 

Here we interpret the X-ray paucity as being due to self absorption. The simplest version of this 
interpretation is to identify the typical GRB peak at ~ 300KeV as being the Wien peak of a thermal 
spectrum with a temperature of ~ 10 2 /T KeV, where T is the bulk Lorentz factor. If we interpret the 
paucity of X-rays from GRB's as being due to a black-body limit, then we can infer that the temperature of 
the radiation as seen by the observer must typically be ~ 10 2 KeV or more. If the brightness temperature 
at the frequency of maximum brightness temperature u m is indeed at the black body limit T m ~ hu m , then 
the energy density in the frame of the fireball, aT^ is limited by the condition that 



A-KR 2 cT 2 aT i < 10 51 L^ergs' 1 . 



(1) 
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The condition that the spectra typically peak at e 7 ~ 300e3oo KeV implies that 

TT m ~ 10 9 ^e 300 . (2) 

This implies that 

r > io 4 - 25 L 5 7i? 13e 2 00 . (3) 

Although this admits unlimited values for both T and R, we note the following simple argument that 
could guide a choice of both: Liberating <~ 10 51 ergs within neutron star dimensions creates fewer photons, 
with larger average energy than dictated by GRB observations, so some additional photon production must 
have taken place. If the (photon creating) dissipation takes place at modest T (say F < 10), then equation 
(3) implies that R < 10 10 cm for L 5 i and e3oo of unity. If the dissipation were to take place at large T and 
perhaps larger R, the question would arise as to why the combination of the two always conspired to yield 
average photon energies of order several hundred KeV. On the other hand, assuming that the dissipation 
takes place at modest T provides a natural reason why the photon production should stop once the average 
photon energy is comfortably below (but not much further below) the pair production threshold: the cycle 
of pair production by photons and photon production by pairs naturally terminates at this point (e.g. 
Cavallo & Rees, 1978; Blandford & Levinson, 1995). For example, we might associate the additional photon 
production with the sharp, large angle deflections the flow might experience while being collimated into jets 
by a surrounding baryonic outflow. If thermalized at modest T, say at a dissipation radius Rdis, and the 
flow expands more or less adiabatically beyond that point, then the temperature would from there on scale 
as 1/R, while T would scale as R as long as the fireball remained baryon pure (to avoid a minimum rest 
mass per unit energy release) and sufficiently optically thick (so that pairs were in thermal equilibrium with 
radiation). Thus equation (3) would continue to be satisfied as long as there were no further processing of 
the thermal photons. Neglecting, tentatively, non-thermal pair production, the thermal pair density would 
become negligible at a temperature well below 3 x 10 8 K (Paczynski 1986; Goodman 1986), and would thus 
establish a photosphere at a radius of about 10Rdi S , with a T of order 10 or less. The value of T could 
increase to 10 3 within a radius of 10 12 cm as long as there were no baryon loading. 

The possibility of non-thermal or grey body (rather than black body) pair production would allow 
for somewhat larger dissipation radius than given by equation (3). This possiblity should not alter many 
of the conclusions of the paper much (though see immediately below) . The basic point is that the modest 
T that we invoke for the production site of the GRB photons seems much smaller than the asymptotic 
values (> 300) that are deduced from other considerations, and the implication is that much of the photon 
production in GRB's comes prior to the acceleration of the flow from T = 10 to F = 300 

Given equation (3), with V of order 10, the characteristic timescale of variation of photons from such a 
photosphere, if of order Rdis/cT 2 , comes out naturally to be of order 3 milliseconds. 

As has been noted many times, gamma rays observed by COMPTEL at about 20 MeV would pair 
produce at these radii unless the bulk Lorentz factor T exceeded about 20. On the other hand, beyond 
the photosphere, the bulk Lorentz factor of a baryon poor wind can grow in proportion to R, so that the 
condition T > 20 is attained not much beyond the photosphere. We contend that any COMPTEL gamma 
ray emitted at R > 3 x 10 10 cm at T > 30 would not be subject to pair production and could in principle 
reproduce time variations as short as 10 ms. 
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In the particularly well studied burst GRB 930131, COMPTEL detected gamma rays out to 20 MeV 
or so with a time profile very similar to that of the BATSE profile, and we suggest that this be interpreted 
as both energy ranges reflecting the intrinsic time profile of the energy release. On the other hand, the 
EGRET gamma rays at 30 MeV, of which there were several, typically arrived several seconds after the 
BATSE and COMPTEL peaks, and this is consistent with the hypothesis that the shortest time profile for 
the EGRET gamma rays was of order seconds. 

3. High Energy Gamma Rays: Production 

Although it is tempting if only for simplicity to model GRB as coming from a particular characteristic 
radius, the following argument suggests that a broad range of emission radii is no less natural. As long as 
the fireball is not baryon loaded and not complicated by non-spherical expansion, the bulk Lorentz factor 
increases linearly with R. The proper time in the frame of the fireball is then described by 

d£ = dt/T oc dR/R. (4) 

It follows that the proper time evolves only logarithmically with R. The duration of shock activity 
(e.g. due to very clumpy baryon contamination, fluctuations in T at launch, etc.) in the proper frame At', 
which is likely to be of order the mean time t' around which such activity is centered, thus persists over 
a large range of radii R. That is, the range of radii over which there is shock activity, Ai?, is given by 
AlnR ~ In(j^), where i?o denotes the beginning of the acceleration zone and might be close to Rdis- If 
» 1, then shocks could persist over many decades of R, particularly shocks associated with the launch 
of the fireball, which can be delayed or "dragged out" to R 3> R by the time dilation. Previous discussions 
of internal shocks occuring at some characteristic radius (Levinson & Eichler 1993, Eichler 1994, Rees & 
Meszaros 1994) have perhaps not emphasized this point. (We agree, however, that shocks resulting purely 
from fluctuations in the saturation radius, at which T approaches its asymptotic value, are likely to form 
near the mean value of that radius if the launch conditions are otherwise identical and if the baryonic 
contamination is not too clumpy.) We contend that there is no reason to expect all of the gamma rays from 
GRB's to come from the same radius. We therefore consider the evaluation of gammaspheric radii as a 
function of gamma ray energy. 

In order to emit gamma rays at energy E, the fireball must be able to accelerate the particles to at 
least energy E, and the gamma rays must then be able to escape. 

The first condition is expressed as 

j m (TTU'/m e c = rieB'/j m m e c (5) 

where 7 m is the maximum Lorentz factor in the frame of the fluid that is achievable by the acceleration, U 1 
is the energy density of the radiation field in the frame of the fireball fluid element, B' is the field strength 
in the fluid frame in cgs units and r\ is the acceleration rate in units of the electron gyrofrequency. For any 
reasonable acceleration mechanism, r\ should be less than unity and for shock acceleration a value of order 
1 to 10 percent seems like a reasonable guess. Writing U' as 10 51 L 51 ergs _1 /47ri? 2 cr 2 equation (5) implies 
that 
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lm /T^0.5^R 12 B'H 51 \ (6) 

For L51 = 1, and assuming that the Poynting flux is limited by 10 51 I/5i ergs per Air steradians, the last 
equation reduces to 

7 m r = 6x 10VHi?f 2 . (7) 

Note that if B = TB' scales as 1 /R, then the maximum energy in the frame of the fireball to which electrons 
can be accelerated increases as i? 1 / 2 T 1 / 2 , as seen from eq. (6). This implies that the ability of the fireball 
to accelerate electrons to extremely high energies increases sharply with radius. At small radii, 

R 12 < 2r]-?B'-? (8) 

a shock cannot accelerate electrons to a high enough energy to pair produce via inverse Compton scattering 
of thermal photons,, (i.e. to a Lorentz factor of 7 m > T, which would be needed given that the photons 
near the spectral peak have energies less than m e c 2 /T in the fluid frame). 

At a radius of order 10 12 cm, with T > 40, equation (7) admits TeV gamma rays which could arrive 
within a fraction of a second if they escape freely. 

If non-thermal pairs are made at larger radii than given in equation (8), they could conceivably 
Compton scatter the thermal photons that were made at smaller radii. If the thermal photons constitute 
most of the burst's photons and most (or much of) its energy, further scattering by pairs downstream 
cannot greatly alter the parameters of the thermal component. 



4. High Energy Gamma Rays: Escape 

In order to compute the pair production opacity contributed by the beamed photons and the associated 
gammaspheric radii, we consider a conical beam of emission with an opening semiangle 9 b , and a power law 
spectrum which is taken to be constant with radius: n(e s ) — n ej a ; e m in < e < e m ax where n is expressed 
in terms of the apparent luminosity, L, as 

l 2 1 -^r 2 

(1 Mb) G (c-max , € m i n ) — — - - "^Gi^rnax •> ^rain) i (9) 



VK^m^r^ / iTT z m e c a r z 

where /j, b = cos6> 6 , G(e max , e min ) = [ln(e mQ;c /e min )] _:L for a = 2, and G(e max , e mm ) = (2 - a)/ (e^ - e^~") 
for «^2. The pair production opacity at energy e 7 (measured in units of mc 2 ) can then be expressed as 

Kp (e 7 ) = 2vr /" de s f d^{(l - n)n(e,)a p (J3)}. (10) 

Here \i m ax — maxj/if,; 1 — 2/e s e 7 }, and <rp(/3) is the pair-production cross section, where p is the speed of 
the electron and the positron in the center of momentum frame, and is given by 

(1 — 1 ) = -T • (11) 

(1 - //)e 7 e s 

The threshold energy, obtained for (3 = and fi = /ifc, now reads: ethrs = 2/(1 — ^&)e 7 — 4T 2 /e 7 . 
Substituting eq. (^|) into eq. ( |lCi| ) yields 

-G{ )A{a : f3 m ax)£-\ (12) 



167rm e c 3 r 2 2 2a r 2Q 
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where 



A(a,f3 max ) 



a + lj 



2\a-l 



(3-/3 4 )ln 



1 + /3 
1-/3 



- 2/3(2 -/3 2 ) 



(13) 



The function A can be evaluated numerically. It vanishes for energies below the threshold energy at which 
Pmax — 0. A plot of A versus a for f3 max — 1 (which is a good approximation well above the threshold 
energy) is given in Blandford & Levinson (1995). Numerically we find A = 0.2 for a photon index a = 2. 



Now, the gamma-spheric radius, r 7 (e 7 ), is defined implicitly by the equation 

/•OO 

/ drn p (r, e 7 ) = 1. 

J r- f 

Using the above expression for the opacity, one obtains 



(14) 



r 7 (e 7 ) = 1.6 x 10 21 L 51 (2T)- 2a G(e max , e mm )A(a, (3 max )e { ^ cm, (15) 

up to a numerical factor that depends on the radial profile of the Lorentz factor T. Adopting, for illustration, 
a = 2, e min = 2, e max = 2 x 10 3 , yields a threshold energy of 10T 2 MeV, a gammaspheric radius 

r 7 (e 7 ) = 2.8 x 10 10 ^-e 7 cm (16) 

above the threshold energy, and a corresponding variability time 

5<= 2^ = 5xl(r5 rf ^ s - (17) 

These last equations admit rather high energy gamma rays emerging from rather compact regions. 
Even the GeV gamma rays observed by EGRET can in principle emerge from within 10 12 cm if the bulk 
Lorentz factor is 300 or more. 

The above results may be considerably altered in the presence of additional radiation component. For 
instance, large angle scattering of the beamed photons by ambient gas surrounding the fireball, e.g., a 
confining, baryon-rich outflow (Nakamura 1998; MacFadyen & Woosley 1998; Eichler & Levinson 1999), can 
produce an unbeamed component that will dominate the pair production opacity within the baryon-poor 
jet even if its luminosity is only a small fraction of the total. It may then be that high energy gamma rays 
can escape only from parts of the baryon-poor jet that have run ahead of surrounding baryons. 



5. Comparison of Constraints on Baryon Contamination 

The cooling of the plasma to 2 x 10 8 K, so that the Boltzman factor for pairs becomes miniscule, is a 
necessary condition for small photospheres, but it is not a sufficient one. Electrons may be entrained from 
the sides (e.g. via neutron seepage) and non-thermal pairs may be produced under some conditions. Let us 
then consider the limits on the net baryon number within a jet that features a compact photosphere. 

Writing 



M = AirmpTiR 2 (3c 



(18) 
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and using cttuR = 1T 2 , where r is optical depth, and letting t be unity at the photosphcric radius R p h, we 
easily obtain 

T a T% h > 2.5 x lO 5 ^!^ (19) 

where T a is defined as L/Mc 2 . 

If the photosphere were to occur at T p h ~ 10, and R p u of order 10 10 cm, then r a would have to be 
enormous, of order 10 6 , and might be convincing evidence that an event horizon enforces baryon purity on 
field lines that thread it. On theother hand, the compactness constraint that motivates our model does 
not imply that the photons that have been produced within 10 10 cm necessarily make their last scattering 
there, it merely implies that they are not swamped by a much larger number of photons made at larger 
radii, nor drained of their energy by mass loading (which would terminate the linear increase of L with R). 
They may be scattered by material that is both dynamically and thermodynamically passive; the timescale 
and spectrum could remain the same. 

We may combine equations (3), (19) and the obvious condition that T a > T p h, to obtain that at 
R P h ~ 10 11 cm, T a must exceed 200. If the photosphere were constrained by some other consideration to be 
larger than this, then equations (3) and (19) similarly constrain T p h to be greater as well. 

6. Conclusions 

We have suggested that a significant fraction of the photons from gamma ray bursts are emitted at a 
radius of about 10 10 cm or less , where the bulk Lorentz factor is modest. This allows the burst entropy to 
be established at modest T, so that the peak energy in the observer's frame can be tied directly to m e c 2 as 
oppsed to a small fraction of m e c 2 to be multiplied by a large T. It also provides a natural account of why 
GRBs have time substructure of milliseconds. 

The photosphere (surface of last scattering) may be larger than the radius of emission, but whether 
inequality (3) is satisfied is independent of R as long as T remains proportional to R. Thus, the fulfillment 
of (3) is unaffected by passive scattering material. 

Higher energy, non-thermal gamma rays may be produced even beyond the photosphere of lower energy 
ones, but we found in section IV that even GeV gamma rays could escape from within 10 12 cm for T of 
order 300. 

TeV gamma rays could, in principle, escape if produced within 10 12 to 10 13 cm for T of order 10 3 
and would be prompt. Thus, it seems worthwhile to run experiments that could detect prompt, ultrahigh, 
non-thermal gamma rays from an unannounced direction. The MILAGRO project is such a detector, though 
the universe may be opaque to the photons (e 7 > 1 TeV) to which it is sensitive. A similar experiment with 
a lower energy threshold - say 10 2 GeV - might detect prompt UHE emission from a sizeable fraction of 
GRB's. Detection of prompt gamma rays above 100 GeV might be the best diagnostic of the bulk Lorentz 
factors of GRB fireballs. 

To conclude, the radius of emission and radius of last scattering of thermal gamma rays in GRB's 
may be different, but both may occur well within 10 12 cm. Non-thermal, high energy gamma rays may be 
produced at and escape from somewhat larger radii than the thermal ones, but they too can emerge from 
compact regions provided that the bulk Lorentz factor is sufficiently high. Soft echos of primary gamma 
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radiation scattered off baryonic matter might provide (via their time scales) limits on the photosphcric 
radius for the bulk of the gamma rays. 

We acknowledge support from the Israel Science Foundation. We thank Dr. E. Waxman for useful 
discussions. 
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